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As part of our comprehensive long-term multi-waveband monitoring of 34 blazars, we followed the activity in 
the jet of the blazar PKS 1510—089 during major outbursts during the first half of 2009. The most revealing 
event was a two-month long outburst that featured a number of 7-ray flares. During the outburst, the position 
angle of optical linear polarization rotated by about 720°, which implies that a single emission feature was 
responsible for all of the flares during the outburst. At the end of the rotation, a new superluminal knot (~ 22c) 
passed through the "core" seen on 43 GHz VLBA images at essentially the same time as an extremely sharp, 
high-amplitude 7-ray and optical flare occurred. We associate the entire multi-flare outburst with this knot. 
The ratio of 7-ray to synchrotron integrated flux indicates that some of the 7-ray flares resulted from inverse 
Compton scattering of seed photons outside the ultra-fast spine of the jet. Because many of the flares occurred 
over time scales of days or even hours, there must be a number of sources of IR-optical-UV seed photons — 
probably synchrotron emission — surrounding the spine, perhaps in a slower sheath of the jet. 
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1. Introduction 



There are four primary methods for probing the 
structure and physics of relativistic jets in blazars 
on parsec and sub-parsec scales: VLBI imaging in 
both total and polarized intensity at millimeter wave- 
lengths, variability of the flux at radio through 7- 
ray frequencies, polarization at radio through opti- 
cal wavebands, and the spectral energy distribution 
(SED). In the past, efforts to study this capricious 
class of objects have been limited by datasets with 
substantial gaps in either time or frequency cover- 
age. Starting in the mid-1990s, this situation has 
been greatly alleviated by the availability of instru- 
ments such as the Very Long Baseline Array (VLBA) , 
the Rossi X-ray Timing Explorer (RXTE), and, most 
recently, the Fermi Gamma-ray Space Telescope. 
Together with concerted efforts with optical/near- 
infrarcd, millimeter-wave, and radio telescopes, these 
facilities now allow long-term, comprehensive multi- 
waveband monitoring of a number of blazars. Such 
programs are providing valuable insights into the pro- 
cesses by which jets form and propagate, as well 
as the locations and physics of flux outbursts in 



blazars 



|e.g 



120081 iMarscher et al. ibooif 



Chatteriee et al.l 120081 lLarionov et al.l 



We are leading a monitoring program that pro- 
vides the data needed to combine the techniques listed 
above in order to locate the sites of high-energy emis- 
sion, determine the processes by which X-rays and 
7-rays are produced, and infer the physical condi- 
tions of the jet, including the geometry of its magnetic 
field. Because of limitations in sensitivity and avail- 
able observing time, such comprehensive monitoring 
can be mounted only for a relatively small number 
of objects, 34 in the case of our program. Progress 
toward an overall understanding of blazars therefore 
involves generalization of inferences drawn from the 
well-studied specimens into a framework that can be 
applied to the entire class. 

We present here a rich set of observations of 
the quasar PKS 1510-089 (z = 0.361) that al- 
lows us to locate the sites of 7-ray flares, leading 
to inferences on the high-energy emission processes. 
PKS 1510—089 possesses all of the characteristics 
of blazars: fiat radio spectrum, apparent superlumi- 
nal motion — among the fastest [as high as 45c for 
H =71 km s~ 1 Mpc~ 1 a nd the concordance cosmol- 
ogy ISj3ergel_^tjiL_ 120071 of all blazars observed thus 
far [Jorstad et al. Il2005f — high-amplitude and rapid 
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flux variability at all wavebands, strong and variable 
optical linear polarization, and high 7-ray apparent 
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luminosity [Hartman et al. Ill999j |. The observations 
of this blazar from our program allow us to compare 
motions, linear polarization, and changing flux of fea- 
tures in the parsec-scale radio jet with flux variabil- 
ity of the entire source at radio, millimeter- wave, IR, 
optical, X-ray, and 7-ray frequencies, and with opti- 
cal polarization. The relative timing of the correlated 
variations thus found probe the structure and physics 
of the innermost jet regions where the flow is acceler- 
ated and collimated, and where the emitting electrons 
are energized. 



2. Observational Results 

Figure[T]presents the Fermi Large Area Telescope 7- 
ray light curve, the R-band optical light curve, and the 
optical polarization as a function of time during the 
first half of 2009, with the data taken from the above 
paper. Detai ls of our observa t ions o f PKS 1510—089 
are given m iMarscher et"al~l |2010} . which includes 
data from 2006 to 2009. That paper also presents 
43 GHz VLBA images showing that a bright superlu- 
minal (22c) knot (A) passed through the core on JD 
2454959±4. Additional VLBA images (see Fig. [2J 
reveal that another knot (B) with the same apparent 
speed was coincident with the core on JD 2455008±8. 
We stress that the "core" of a blazar seen on VLBI 
images lies parsecs downstream of the location of the 
black hole [Marscher et al. 1 |2002[ | . There is mount- 
ing evidence that, barring optical depth effects, the 
core represents emission from jet plasma that has been 
compressed and energized by a system of one or more 
standing conical shocks owing to pressure mismatches 
between the jet and the external medium. 

One of the main questions that our observational 
program can answer is whether the high-energy emis- 
sion occurs between the central engine and core, 
within the core, or downstream of the core. In or- 
der to do this, we must first establish an association 
between a superluminal knot and high-energy flares, 
then analyze the timing of the flares relative to the 
motion of the knot. In the case of BL Lac in 2005, 
we employed this method to determine that there are 
multiple sites of X-ray and optical flares in the jet: 
upstream of the core where the magnetic field has a 
helical geometry, and within the core. Flares occur as 
a disturbance in the jet flow (an "emission feature" ) 
passes through these regions [Marscher et al. |[2008j [ . 

Eight major flares are apparent in Figure [TJ which 
covers the first half of 2009. The figure numbers the 
flares. (We note that 7-ray emission was detected at 
an elevated level throughout the 50+ days following 
flare 8, but one cannot isolate distinct flares during 
this period.) Of particular interest is the 51-day pe- 
riod from JD 2454901 to 2454962 during which flares 
3-8 took place. As seen in Figure [TJ the optical po- 
larization position angle x rotated by ~ 720°. (Since 



there is no distinction between \ an d X ± 180°n, we 
select n such that the jump in x is minimized.) The ra- 
tio of 7-ray 0.1-200 G eV integrated flux to b olometric 
synchrotron flux [see IMarscher et al. Il2010l | for flares 
1-8 was 70, 30, 40, 40, 30, 10, 40, and 9, respectively. 

Superluminal knot A passed through the core at the 
same time (within the uncertainties) as the very sharp 
(intra-day), extremely high- amplitude 7-ray and op- 
tical flare on JD 2454962. The rotation of the opti- 
cal polarization vector ended at the same time. We 
therefore surmise that the 7-ray flares occurred up- 
stream of and within the core. After this point, the 
7-ray flux remained elevated, which suggests that 7- 
ray emission continued downstream of the core, as in- 
ferred previously from observations wit h EGRET, the 
VLBA, and single-dish radio antennas [Jorstad et al. I 
l200li lLatheenmaki fc Valtaoia 112003] . 



3. Interpretation 

Although a turbulent magnetic field can produce 
such apparent rotations of the pola rization vector by 
720° or more via a ran dom walk lAller et al.l Il985l 
iD'Arcangelo et all 120071 iJonesI [1988|, the beginning 
and end of such a rotation should occur randomly, 
whereas the observed event matches quite closely the 
duration of the flux outburst. We therefore conclude 
that a single coherent emission feature was responsible 
for the entire ou tburst. We appeal to t he same model 
as for BL Lac [Marscher et al. 1 12008} : as a distur- 
bance in the flow moves down the jet, it follows a spiral 
streamline [as predicted for a model in whi ch the jet 
is ma gnetically accelerated and collimated; IVlahakisl 
2006], It first travels through the acceleration and col- 
limation zone, which possesses a tightly wound helical 
magnetic field. If the emission feature covers most, 
but not all, of the cross-section of the jet, much of 
the polarization is cancelled out from mutually per- 
pendicular magnetic field orientations around the cir- 
cumference. The residual polarization vector rotates 
as the feature executes its spiral motion. 

We have performed a rough calculation of the rota- 
tion of x m such a model, matching the parameters 
(probably non-uniquely) to the PKS 1510—089 data. 
The bulk Lorentz factor increases from 8 to 24 dur- 
ing the rotation event, causing the Doppler factor to 
increase from 15 to 38 for an angle between the je t 
axis and the line of sight of 1.4° [Jorstad et al. 2005]. 
The acceleration of the jet causes an increase with 
time in the rate of rotation, as observed. The circula- 
tion of the centroid of the emission feature around the 
jet axis causes x (corrected for relativistic aberration) 
to undulate with time about the trend of monotonic 
rotation. In the quantitative model, the core lies 17 
pc downstream of the point where the outburst con- 
taining flares 3-8 began. The magnetic field is in the 
range of 1 to 0.2 G during the outburst. 
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The dichotomy in the ratio of 7-ray to synchrotron 
flux of the various flares suggests that different sources 
of seed photons dominate the inverse Compton scat- 
tering as the emission feature moves down the jet. 
Flares with relatively low values of this ratio can be 
caused by an increase in the number of GeV electrons, 
which causes enhanced synchrotron, synchrotron self- 
Compton (SSC), and external Compton (EC) emis- 
sion. If SSC is the main mechanism, second-order 
scattering is probably important, since the lowest 7- 
ray to synchrotron flux ratio, 9, still exceeds unity. 
There are some flares (1, 3, and 7) where there is lit- 
tle or no increase in the optical flux. We infer that 
these correspond to the emission feature encountering 
local sources of enhanced IR-optical-UV seed photon 
density. We estimate the luminosity of these sources 
to be ~ 3 x 10 43 erg/s, too high to correspond to any 
likely cosmic source except a slower (but still proba- 
bly r elativistic) sheath of the jet [see iGhisellini et al.l 
2005]. The source could, for example, be a shock in 
the outer periphery of the jet. 

Figure [3] sketches the general picture that we pro- 
pose. The 7-ray emission undergoes flares at various 
sites along the jet, with events 3-7 representing loca- 
tions where either the sheath produces large numbers 
of seed photons or electrons are accelerated to GeV 
energies. Flare 8 takes place as the moving emission 
feature becomes compressed and energized by stand- 
ing shocks in the core. The elevated post-outburst 7- 
ray emission first comes from knot A as it propagates 
downstream of the core. At some point — probably 
near JD 2454990 — emission from knot B becomes 
dominant. This is weaker than knot A, so the flares 
are of lower amplitude. But we do see a ~ 180° 
rotation of the polarization vector in the same direc- 
tion and with the same slope as the previous, longer 
rotation. This repetition of the pattern is a firm re- 
quirement of our model: the physical structure of the 
spine or sheath of the jet could change on a time scale 
of years, but not months. In contrast, a turbulent 
field model predicts that both the sense and rate of 
rotation should vary from one event to another. 

It is possible that flares 1-2 correspond to the distur- 
bance passing through the rich seed photon field near 
the accretion disk and within the broad emission-line 
region (BELR). If this is the case, then there is an 
extended region in the jet beyond the BELR where 
electrons are not accelerated efficiently to energies ex- 
ceeding ~ 1 GeV. 

4. Conclusions 

We are finally at the point where the richness of our 
datasets is sufficient for us to draw grand inferences 
about the locations and physical mechanisms of 7-ray 
flares. If PKS 1510-089 and BL Lac are typical, it 



is no wonder that previous, less comprehensive moni- 
toring programs left us confused! There are multiple 
sources of seed photons, some of which are quite local 
(as opposed to more global, such as an IR-emitting 
hot dust torus) and may not radiate at sufficiently 
high luminosities to be directly observable. These in- 
clude the accretion disk and BELR — which may be 
important during the earliest part of a multi-flare out- 
burst — shocks or other features in the sheath of the 
jet, and synchrotron radiation from the fast spine of 
the jet that contains the electrons that scatter the seed 
photons to 7-ray energies. A single superluminal knot 
can be responsible for a number of flares and periods 
of sustained elevated 7-ray emission. 

As we accumulate data for more blazars, we expect 
to see similar patterns of flares both before and af- 
ter a disturbance passes through the core of the jet. 
We also hope to find deviations from this pattern that 
serve to provide further insight into the range of phys- 
ical behavior of blazars. In both cases, we antici- 
pate a tremendous increase in our understanding of 
the relativistic jets of blazars by combining the great 
power of monitoring the flux, polarization, and sub- 
milliarcsecond structure of these exciting objects. 
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Figure 1: Variation with time during the first half of 2009 of (top to bottom panels) (1) 7-ray flux (upper limits in 
red), (2) optical flux density (note the logarithmic scale, to fit in the very high-amplitude flare on JD 2454962), (3) 
degree of optical linear polarization, and (4) electric-vec tor position angle of op tical linear polarization. Vertical dotted 
lines denote major 7-ray flares. From data presented bv lMarscher et al. I |2010l ]. 
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Figure 2: VLBA image of PKS 1510-089 on 16 September 2009 (JD 2455091). The resolution beam is comparable to 
the angular size of the longest baselines of the array in the direction of the jet. Contour levels are -1, 1, 2, 4, 8, 16, 32, 
and 64% of the peak intensity of 0.55 Jy/beam. The two superluminal knots discussed in the text are marked, as is the 
core. 
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Figure 3: Sketch of the model for PKS 1510—089 discussed in the text. Note the logarithmic length scale, so that 
various components of the nucleus can be included. 
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